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SECTION  i 


INTRODUCTION 


A gas  dynamic  'laser  (GDL)  Is  a combustoi'-driven  device  that  requires 
propellants  for  operation.  The  propellants;  are  supplied  to  the  combustor  by 
the  fluid  supply  system  (FSS). 

The  FSS  is  a series  of  tanks  (three  propellant  tanks  and  two  prrssurant 
tanks),  associated  plumbing,  and  a control  system  to  regulate  the  propellant 
flow  and  tank  pressures.  The  propellant  tanks  are  pressurized  with  a helium 
pressurant.  The  helium  pressurant  supplies  the  energy  for  the  propellant  trans- 
port from  the  tanks  to  the  combustion  chamber. 

The  r jlium  pressurant  tanks  are  pressurized  from  an  ambient,  high-pressure 
helium  source  tank.  The  helium  from  the  source  tank  Is  pressure  controlled  to 
the  lower  pressure  of  the  pressurant  tank  and  cooled  to  near  liquid  nitrogen 
temperatures  by  a nitrogen  bath  heat  exchanger.  The  pressurant  tanks  are  double- 
wall vacuum  jacketed  tanks  with  external  cooling  coils. 
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SECTION  II 

ANALYSIS  OF  THE  FSS  LOADING  SYSTEM 

A schematic  of  the  pressurant  loading  system  Is  shown  In  figure  1. 

It  can  be  shown  that  the  temperature  of  the  pressurant  helium  In  the  FSS  at 
the  start  of  the  operation  of  a GDL  should  be  nearly  the  same  temperature  as  the 
propellant  It  will  pressurize  for  best  operation.  This  procedure  reduces  the 
heat  and  mass  transfer  between  the  pressurant  and  the  propellant  In  the  propel- 
lant tank.  This  pressurant  temperature  condition  requires  one  pressurant  tank 
to  be  cryogenic  (140°R)  and  the  second  to  be  near  ambient  (400°R). 

This  report  describes  the  analysis  of  the  pressurant  loading  system  to 
determine  the  temperature-,  pressure-,  and  mass-time  history  of  the  helium 
pressurant  system. 

The  pressurant  loading  of  the  two  tanks  can  be  conducted  several  ways. 

These  include:  (1)  simultaneous  tank  loading  (load  both  tanks  at  the  same  time), 

(2)  sequential  loading  (load  one  tank,  then  transfer  its  contents  to  the  second 
tank  while  maintaining  the  first  tank  at  constant  pressure,  and  (3)  Individual 
loading  (load  each  tank  independently). 

The  analysis  of  the  pressurant  loading  systems  was  made  by  taking  energy  and 
mass  balances  around  each  tank. 

Time  dependent  energy  equation: 

(me)  - 6Q  - 6Wn  + £(hm)j  - £(hm)k  (1) 

j k 

Time  dependent  mass  equation: 

at  <">  * - £ mk  <2> 

Heat  transfer  between  the  helium  and  the  tank  walls  was  considered.  The 
temperature  distribution  within  the  wall  was  determined  from  the  time  dependent 
heat  conduction  equation. 
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Figure  1.  Pressurant  Fill  Process  Schematic 
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VTm  n 32T 

a — 7 


(3) 


The  initial  and  boundary  conditions,  as  well  as  the  additional  equations, 
necessary  for  solution  of  the  set  of  differential  equations  are  determined  from 
the  constraints  Imposed  on  the  system. 

Several  assumptions  have  been  made  to  simplify  the  analysis. 

1.  Instantaneous  and  homogeneous  mixing  of  any  Incoming  helium  streams 
and  the  helium  contents  of  the  artks. 

2.  No  energy  or  pressure  losses  due  to  flow-through  Interconnecting 

1 ines. 

3.  The  pressure  controllers  operate  Isenthalplcly. 

4.  The  heat  exchanger  Is  100  percent  efficient. 

5.  The  cooling  coils  on  the  pressurant  tanks  distribute  their  effect 
over  the  entire  surface  of  the  tank,  not  just  point  locations. 


6.  The  wall  temperature  distribution  can  be  approximated  by  a third 
order  polynomial  in  the  spatial  coordinant. 

The  set  of  differential  equations  was  evaluated  numerically  using  a modified 
Adams  method  as  implemented  in  DIFSUB*  on  the  CDC  6600  at  t.he  Air  Force  Weapons 
Laboratory,  Kirtland  AFB,  New  Mexico.  The  necessary  helium  thermodynamic  data 
was  supplied  by  the  National  Bureau  of  Standards  and  Is  that  found  in  NBS  Tech- 
nical Note  622.** 


* J.  W.  Aubrey,  Supporting  Mathematical  Group,  Technology  Division,  Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico 

**  Robert  D.  McCartey,  Cryogenic  Division,  Institute  of  Basic  Standards, 
National  Bureau  of  Standards,  Boulder,  Colorado 
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SECTION  III 
SAMPLE  RESULTS 

Although  the  set  of  differential  equations  could  be  used  to  analyie  the 
several  different  pressurant  loading  schemes,  the  most  promising  $cha..«e  appears 
to  be  the  sequential  loading.  In  this  procedure,  helium  Is  loaded  Into  one 
tank  with  precnllled  walls.  This  tank  will  be  the  cryogenic  temperature  tank. 
Once  the  pressure  In  the  cryogenic  temperature  tank  reaches  a set  value,  e.g., 
the  maximum  pressure,  helium  Is  allowed  to  flow  from  the  cryogenic  temperature 
tank  to  the  second  (ambient  temperature)  tank.  The  walls  of  the  ambient  tem- 
perature tank  may  or  may  not  be  Initially  at  ambient  temperatures.  Flow  Into  * 
the  ambient  temperature  tank  Is  controlled  so  that  a constant  pressure  Is  main- 
tained within  the  cryogenic  temperature  tank.  During  the  loading  operation, 
the  cooling  coll  on  the  cryogenic  tank  Is  removing  energy  from  the  helium  In 
that  tank.  Once  the  maximum  pressure  within  both  tanks  Is  reached,  the  helium 
flow  between  the  tanks  Is  terminated.  The  cryogenic  temperature  tank  is  cryo- 
pumped  to  Its  final  temperature,  where  upon  the  pressurant  tank  loading  Is 
completed. 

Figure  2 shows  the  tempera ture-tlme  history,  figure  3 shows  the  pressure- 
time history  and  figure  4 shows  the  mass-time  history  for  the  two  pressurant 
tanks  for  a sequential  loading  procedure. 
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TEMPERATURE  VS  TIME 


Figure  2.  Predicted  Temperature  History  of  Pressurant 


PRESSURE  VS  HUE 


Figure  3.  Predicted  Pressure  History  of  Press urant  Tank 


MASS  VS  TIME 


Figure  4.  Mass  History  of  Pressurant  Tanks 
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APPENDI/ 

DEVELOPMENT  OF  EQUATIONS 


(4) 


3I=  ^mj  ' ^mk  <5) 

Equations  (4)  and  (5)  describe  the  energy  and  mass  time  rates  of  change, 
respectively,  for  any  control  volume  1.  For  clarity  the  i subscript  has  been 
omitted.  In  addition,  two  definitions,  equations  (6)  and  (7)  are  presented. 


H = E + PV 

(6) 

> 

CL 

II 

E 

(7) 

Equation  (6)  defines  enthalpy;  equation  (7)  defines  density.  If  only 
internal  energy  is  significant,  equation  (6)  can  be  rewritten. 


mh  = me  + PV 


(8) 


Time  derivatives  can  be  taken  of  equations  (8)  and  (7). 


mdh  + 

Ht 


hdm  _ d 

HT  HT 


(m§)  + 


Pd V 
dt 


VdP 

nr 


(10) 


n 


For  the  analyris,  energy  and  mass  balances  are  taken  on  the  helium  within 
well  defined  control  volumes,  i.e.,  the  tanks.  Since  helium  is  a pure  material 
and  will  be  only  in  the  gaseous  phase,  the  Gibbs  phase  rule  specifies  any  inten- 
sive property  can  be  determined  by  two  other  intensive  properties,  specifically 
temperature  and  pressure.  Specific  enthalpy  and  density  can,  therefore,  be 
expressed  as  functions  of  temperature  and  pressure  only. 


h = h (T,  P) 
P ■ P (T,  P) 


(ID 

(12) 


Substituting  equations  (10)  and  (11)  into  (4),  substituting  equation  (12) 
into  (9),  and  rewriting  equation  5,  yield  the  following  equations. 


m9h 

at 


dT  + 

at + 


(T!t  - v)  Sj|  = «Q  - «Hn  ♦ Etta) . - £(B)k  (13) 


V3p 

aTl 


41+  V3p 
dt  5TH 


dP  pdV  dm  _ n 

dt  dt  ’ “ dt  " 


dt 


m, 


k 


(14) 


(15) 


Equations  (13),  (14),  and  (15)  form  the  basis  for  the  analysis  of  the  helium 
loading  system  for  the  FSS.  In  addition  to  the  three  equations,  four  more 
relations  are  required. 


v<t> 


Q = Q(t) 


wn  = W(t) 


m = m(t) 


12 


(16) 

(17) 

(18) 

(19) 
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These  four  equations  are  obtained  after  determlng  the  constraints  that  are 
assumed  to  apply  to  the  system.  For  the  analysis 


V'(t)  = 0 


W(t)  * 0 


These  conditions  are  due  to  the  rigid  wall  construction  of  the  propellant 
tanks. 

mft)  = arbitrary  constant 

This  condition  is  dependent  upon  the  capacity  of  the  heat  exchanger  and  Is 
a necessary  but  arbitrary  function  for  the  solution  of  the  system  of  equations. 

Of  importance  is  the  heat  rate  equations  for  the  cryogenic  and  ambient  tem- 
perature helium  tanks.  Since  the  tank  walls  are  a large  thermal  mass,  heat  flux 
to  and  from  these  walls  can  greatly  effect  the  helium  conditions  within  the 
tanks. 

Heat  is  assumed  to  flux  through  the  tank  walls  due  to  tank  interior  and 
exterior  temperature  differences.  Convective  or  Newtonian  heat  transfer  is 
assumed  at  the  tank  wall  surfaces.  Conductive  or  fourier  heat  transfer  is 
assumed  within  the  tank  wall. 

The  tanks  of . the  FSS  are  spherical.  Since  the  tank  radius  is  generally  much 
greater  than  the  tank  wall  thickness  and  the  heat  flux  is  assumed  to  be  evenly 
distributed  on  the  tank  wall  surface,  the  heat  conduction  equation  can  be  written 
approximately  in  one  dimensional  cartessian  coordinates. 


The  initial  condition  is 


TJC.O)  = Tw1  alt  5 
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.'»*«*** 


The  boundary  conditions  are 


F*  (L.t)  = Jji  (T-Tw(L,t)) 


t > 0,  £ - L 


jf  (0,t)  = jS.  (Tw(0,t)  - Tc)  t > 0,  £ - 0 


The  Initial  and  boundary  conditions  are  not  homogeneous  so  an  analytic 
solution  by  separation  of  variables  Is  not  possible.  The  wall  temperature  dis- 
tribution can  be  approximated  by  assuming  a polynomial  T = P(£,t)  where  the 
coefficients  are  to  be  determined.  For  the  analysis  a third  order  polynomial 
in  £ was  assumed. 


Tw  r B£3  + 0£2  + D£  + E' 


Substituting  equation  (24)  into  (22)  and  (23)  lead  to  two  algebraic  equa- 
tions. Substituting  equation  (24)  Into  (20)  and  solving  for  £ * 1/3L  and  £ a 2/3L 
leads  to  two  differential  equations  for  the  time  rate  of  change  of  the  coeffi- 
cients B and  C of  the  form 


c M + f 3C  + f 3T+f  °'c  _ r 

i at  2 at  r3  at  4 at  " 


n I C 

is  obtained  from  equations  (13)  through  (15).  ~9y  , the  time  rate  of 
change  in  the  coolant  temperature  has  to  be  approximated. 


0'r  i 2 

r=l0“r  <Tc  - V 


This  arbitrary  approximation  assumes  that  the  response  of  the  coolent  tem- 
perature is  100  times  greater  than  the  temperature  response  of  one  side  of  the 
tank  wall  to  a temperature  change  to  the  other  side. 

* 

Once  the  temperature  distribution  is  known  the  heat  rate  6Q  can  be  evaluated. 


6Q  = njA  (T  - B£3  - C£2  - D£  - E ’ ) 
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.-'urther  the  heat  load  on  the  coolant  coll  can  be  evaluated. 

Qc  B n0A  ( E * - Tc)  (27) 

Tne  heat  transfer  coefficient  n j was  approximated  by  evaluating  the  natural 
convective  heat  transfer  coefficient  for  a verticle  flat  plate  with  a length 
equal  to  the  tank  diameter.  The  fluid  properties  were  evaluated  assuming  the 
pressure  to  be  the  maximum  and  the  temperature  to  be  the  time-averaged  tempera- 
ture for  the  helium  within  the  tank. 

The  heat  transfer  coefficient  nQwas  approximated  b.>  evaluating  the  heat 
transfer  coefficient  for  turbulent  flow  in  tubes5  i.e.,  the  cooling  coll,  for 
a constant  heat  flux.  The  coefficient  obtained  was  then  reduced  by  the  ratio  of 
cooling  coil  shadow  area  to  the  total  tank  area.  The  fluid  properties  were 
evaluated  for  liquid  nitrogen. 
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Variables 


A 

B 

C 

D 

e 

E 

E' 


G 


h 

H 

k 

L 


m 

m 


P 

Q 


T 

T 


c 


t 

V 
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ABBREVIATIONS  AND  SYMBOLS 


Meaning  Units 


surface  area  of  tank  for  heat  flux  ft2 

first  coefficient  for  wall  temperature 

distribution  °R/ft3 

second  coefficient  for  wall  temperature 

distribution  °R/ft2 

third  coefficient  for  wall  temperature 

distribution  °R/ft 

specific  internal  energy  BTU/lb 

energy  BTU 

fourth  coefficient  for  wall  temperature 
distribution  °R 


function  Involving  the  heat  transfer  coefficients 
and  wall  thickness,  used  in  equation  (25)  to 
evaluate  the  wall  temperature  distribution 
coefficients 

function  involving  the  thermal  diffusivity  and 
wall  thickness,  used  in  equation  (25)  to  evalu- 
ate the  wall  temperature  distribution  coefficients 


specific  enthalpy  BTU/lb 

enthalpy  BTU 

thermal  conductivity  BTU/sec-ft-°R 

wall  thickness  ft 

mass  BTU 

mass  rate  lbs/sec 

pressure  psia 

heat  rate  BTU/ sec 

helium  temperature  °R 

coolant  temperature  °R 

wall  temperature  for  equation  (21)  °R 

time  sec 

volume  ft3 

work  rate  BTU/sec 

thermal  diffusivity  ft2/sec 
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ABBREVIATIONS  AND  SYMBOLS  (cont'd) 


Meaninc 


heat  transfer  coefficient  of  tank  interior  wall 
heat  transfer  coefficient  of  tank  exterior  wall 
length  along  one-dimensional  coordinant 
density 


Units 

BTU/sec-ft2-°R 

BTU/sec-ft2-°R 

ft 

lb/ft3 


denotes  initial  conditions 
denotes  flow  into  a control  volume 
denotes  flow  out  of  a control  volume 
denotes  constant  pressure 
denotes  constant  temperature 


summation 

change  in;  AT  = T2-Ti 
time  derivative 

partial  derivative  with  respect  to  pressure 

partial  derivative  with  respect  to  temperature 

partial  derivative  with  respect  to  time 

partial  derivative  with  respect  to  one-dimensional 
coordinate  £ 

second  partial  derivative  with  respect  to  one- 
dimensional coordinate  5 


l 


